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(54) Tide: VEHICLE OCCUPANT DISCRIMINATION SYSTEM AND METHOD 




(57) Abstract 

An occupant discrimination system incorporates a position sensor (16) for measuring the distance from a fixed structure (28) within 
the vehicle (10) to a first surface (120) of an occupant (24) or object on a vehicle seat, and a seat weight sensor (352) for measuring the 
weight of the occupant or object The position sensor (16) generates a plurality of signal components that are combined with the weight 
and vehicle speed related measurements so as to form a signal space, from which a plurality of measurements are calculated. A plurality of 
seat occupancy scenarios are defined and the likelihood of each seat occupancy scenario is functionally related to the plurality of measures. 
For a given measure state, the likelihood of each of the seat occupancy scenarios is calculated, and the most likely scenario is used to 
govern the control of a safety restraint system (38). 
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VEHICLE OCCUPANT DISCRIMINATION SYSTEM AND METHOD 



TECHNICAL ART 



The instant invention generally relates to systems and methods for discriminating the 
seat occupancy configuration in a vehicle for purposes of controlling the activation of a 
5 safety restraint system. 



BACKGRO UND OF THF. INVENTION 

Recently, the automotive industry has been faced with the challenge of designing a 
vehicle airbag deployment system which can inhibit deployment of the airbag in situations 
where allowing deployment increases the likelihood of causing harm to an occupant, such as 
where a rear-facing infant seat is located in the front passenger seat. Examples of such 
attempts or proposals include: (1) manually actuated switches for enabling/disabling air bags; 
(2) magnetic or inductive "tags" that must be affixed to a rear-facing infant seat; and (3) 
capacitance sensors to distinguish humans from inanimate objects, or weight threshold 
sensors that disable the air bag below a predetermined threshold weight (such as 25 pounds). 

However, none of these arrangements has proven entirely satisfactory. More 
specifically, manually actuated switches are subject to human forgetfulness or inattention. 
Tag systems have to be retroactively affixed to all existing rear-facing infant seats. 
Capacitive and weight sensing systems have not achieved acceptable levels of reliability. 
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SUMMARY OF THE INVENTION 

The instant invention overcomes the above-noted problems by providing a vehicle 
occupant discrimination system and method which is capable of differentiating rear-facing 
infant seats from people so that deployment of passenger-side air bags can be controlled 
5 depending on whether a rear-facing infant seat or larger human is present. 

The system and method of the instant invention reliably distinguishes humans from 
infant seats without any action on the part of the driver or other person. The instant 
invention examines a variety of information, such as from separate optical, weight, speed, 
and door switch sensors, to determine factors such as distance between the airbag and the 
10 object located in the seat, small amplitude motions such as breathing, large scale motions, 
and weight. The instant invention tracks current and past information concerning these 
factors to accurately discriminate between people and infant seats. 

In accordance with the instant invention, a vehicle occupant discrimination system and 
method comprises an algorithm programmed into processor that processes information 

15 produced by an optical sensor, such as a sensor utilizing a collimated and modulated 
optical beam produced by an LED, which is reflected from a surface and imaged by a lens 
onto a position sensitive detector (PSD). A range limit detector can also be utilized to 
provide an indication when the reflecting surface is too close to the PSD to produce a 
reading. The position of the imaged light on the PSD is used to compute distance to the 

20 reflecting surface from a fixed structure within the interior of the vehicle. As discussed 
more fully hereinbelow, the instant invention also preferably processes information from 
other inputs such as a seat weight sensor, a speed sensor, and a door switch. 

The system and method of the instant invention tracks distance, patterns of motion, and 
weight over time and processes to determine the likelihood of various predefined seat 
25 occupancy scenarios. As a particular scenario becomes "older," the discrimination thereof 
becomes more reliable. Examples of particular seat occupancy scenarios include an empty 
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seat, a seat with a rear facing infant seat (RFIS), a normally seated occupant, an a normally 
seated occupant holding a newspaper or map. 

The signals from the optical sensor are identified as follows: 

1 . Inside channel signal (I): this is the signal from the inside end of the PSD chip. 
The signal level gets higher as the distance to the target gets higher; 

2. Outside channel signal (O): this is the signal form the outside end of the PSD 
chip. The signal level gets lower as the distance to the target gets higher; 

3. LED signal: this is the value of the driving voltage of the LED; and 

4. Limit detector signal: this switch signal becomes high when the target is closer 
than a threshold. 

In addition, processor receives weight information from the weight sensor and speed 
information from the speed sensor. The door switch is used for instantaneous recognition 
of new scenarios. 



In accordance with a first and preferred embodiment of the instant invention, the 
15 following four measures are calculated from the input signals: 

(M 0 ) Mean distance (avgda); 
(M,) Variance of the distance (a 2 ); 
(M 2 ) Mean weight (w); and 
(M 3 ) Slant measure (8). 

20 These measures are reset at the onset of a new scenario and therefore the recognition of 

a new scenario is a key factor. 
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An associated scenario sample space {S,} is divided into the following three categories: 

(S 0 ) Rear Facing Infant Seat (RF1S); 
(S,) Empty Seat; 

(S 2 ) Person Seated Normally (PSN); 

5 In accordance with a second embodiment of the instant invention, the following five 

measures are calculated from the input signals: 

(M 0 ) Mean distance (avgda); 
(M,) Maximum distance (Xmax); 
(M 2 ) Variance of the distance (a 2 ); 
10 (M 3 ) Mean weight (w); and 

(M 4 ) Slant measure (6). 

The assoicated scenario sample space {SJ is divided into the following eight 
categories: 

(S 0 ) Rear Facing Infant Seat (RFIS); 
15 (S,) Front Facing Infant Seat (FFIS); 

(5 2 ) Object; 

(5 3 ) Empty Seat; 

(5 4 ) Rear Facing Infant Seat with Blanket (RFIS/b); 

(5 5 ) Person Seated Normally (PSN); 

20 (S 6 ) Person Reading Newspaper (PRN); and 

(S 7 ) Person Holding Object (PHO); 
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All of these scenarios are assigned a priori probabilities Pr^). In the case of the second 
embodiment, the probability distributions defined for the above five measures given these 
eight scenario, i.e., f x (M$ S„ thus making a total of 40 distributions. As processor 
calculates the five measures, probability of a particular measure given a particular scenario 
5 (i.e., Pr(Mj Sj)) is calculated using the available distributions. The probability of the five 
measures, given a particular scenario, is calculated as follows: 

4 3 4 

Pr(M 0l M ll M 2 .M 3 .M 4 |S i ) = n p KMj|Si)Wji £ Pr(Mj|Si)Pr(Mk|Si)Pr(Si){C V1C ) 

j*k, k>j 

10 where, W j{ and C jkKi are the autocorrelation and crosscorcelation terms respectively. The 

correlation matrix (Ij ) for each scenario for the five measures is : 
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It is assumed that there is very little correlation between the measures, therefore the 
crosscorrelation terms are taken to be zero. The autocorrelation terms being one, the above 
equation then simplifies down to : 



Pr(M 0t M 1 ,M 2 ,M 3l M 4 |S i ) = r7Pr(Mj|Si) 

j-0 

Then, probability of a scenario can be calculated using the five measures and Baye's 
20 rule as follows: 



PrrciM m m m ■ Pr(M o,M i,M 2,M 3,M 4lSi)Pr(SQ 

PrtSjlMo.M^Mj.Ma.M^ = —7 - - - 

X Pr(M o,M i,M 2 ,M 3,M 4|Sk)Pr(Sk) 



k-0 
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This provides eight probabilities, ie, probabilities of the eight scenarios, given the five 
measures, Pr(S i |M 0 ,M J ,M 2 ,M 3 ,M4,M 5 ) e. The airbag inflator is then enabled, controlled, or 
disabled depending on what the most probable scenario warrants. 

Accordingly, one object of the instant invention is to provide an improved vehicle 
5 occupant discrimination system that can distinguish various seat occupancy scenerios. 

A further object of the instant invention is to provide an improved vehicle occupant 
discrimination system that can be used to affect the control of a safety restraint system so 
as to reduce the risk of injury to occupants in the event of a crash or the deployment of the 
safety restraint system. 

10 These and other objects, features, and advantages of the instant invention will be more 

fully understood after reading the following detailed description of the preferred 
embodiment with reference to the accompanying drawings and viewed in accordance with 
the appended claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 FIG. 1 is a schematic illustration of vehicle occupant discrimination system having an 

occupant position detector and a weight sensor in accordance with the present invention 

FIGS. 2a, 2b and 2c are schematic illustrations of an embodiment for detecting the 
occupant position. 

Fig. 3 illustrates a block diagram of the hardware incorporated in various embodiments 
20 of the instant invention. 

Figs. 4a, and 4b illustrates an overall block diagram of a signal processing algorithm in 
accordance with a first embodiment of the instant invention. 
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Fig. 5 illustrates a block diagram of a subprocess for detecting if the signal from the 
position sensor has a low level, in accordance with the process of Figs. 4a, and 4b. 

Fig. 6 illustrates a block diagram of a subprocess for calculating the distance to an 
occupant or object on the seat of a vehicle, in accordance with the process of Figs. 4a, and 
5 4b. 

Fig. 7 illustrates a block diagram of a subprocess for calculating the slant component in 
the measurement of distance to the occupant or object on the seat of a vehicle, in 
accordance with the process of Figs. 4a, and 4b. 

Fig. 8 illustrates a block diagram of a subprocess for filtering vehicle speed and 
10 calculating and filtering acceleration from the vehicle speed signal, in accordance with the 
process of Figs. 4a, and 4b. 

Fig. 9 illustrates a block diagram of a subprocess for detecting big changes in the 
distance measurement, in accordance with the process of Figs. 4a, and 4b. 

Fig. 10 illustrates a block diagram of a subprocess for calculating primary 
15 measurements, in accordance with the process of Figs. 4a, and 4b. 

Fig. 11 illustrates a block diagram of a subprocess for calculating secondary 
measurements, in accordance with the process of Figs. 4a, and 4b. 

Fig. 12 illustrates a block diagram of a subprocess for handling low level signal 
conditions, in accordance with the process of Figs. 4a, and 4b. 

20 Fig. 13 illustrates a block diagram of a subprocess to check if the seat occupancy 

scenario is settled, in accordance with the process of Figs. 4a, and 4b. 

Fig. 14a, 14A, and 14c illustrates a block diagram of a subprocess for making primary 
probability calculations, in accordance with the process of Figs. 4a, and 4b. 
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Fig. 15a, ISA, and 15c illustrates a block diagram of a subprocess for making 
secondary probability calculations, in accordance with the process of Figs. 4a, and 4b. 

Fig. 16a, 16b, and 16c illustrates a block diagram of a subprocess for checking for 
resets to the seat occupancy scenario, in accordance with the process of Figs. 4a, and 4b. 

5 Fig. 17 illustrates an overall block diagram of a signal processing algorithm in 

accordance with a second embodiment of the instant invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Referring to FIG. 1, a vehicle occupant discrimination system 10 in accordance with 
the instant invention senses and discriminates the condition of seat occupancy by an object 

10 on a vehicle seat for controlling the activation of a safety restraint system, such as an air 
bag inflator 38 or for activating a warning device 40. More particularly, the vehicle 
occupant discrimination system 10 comprises a position sensor 300 and a seat weight 
sensor 400. The position sensor 300 measures the distance of an object in the vehicle seat 
relative to the fixed structure within the vehicle and generates transitory values for 

15 occupant position therefrom. The weight sensor 400 generates a weight output 406 
responsive to the weight of an occupant 24 or object on the vehicle seat 28. 

In one embodiment of the instant invention as illustrated in Fig. 1, the weight sensor 
400 comprises a fluid filled bladder 402 together with a pressure transducer 404 for 
measuring the pressure therewithin. However, the instant invention is not limited to any 
20 particular design or configuration of the associated weight sensor 400 in the seat so long as 
the particular weight sensor 400 provides a measure of the weight of the occupant 24 or 
object on the vehicle seat 28. 

In one embodiment of the instant invention, the position sensor 300 comprises active 
infrared position detector 12 comprising the following : a modulator 14 connected to an 
25 optical transmitter 16, and an optical receiver 18 connected to a synchronous 
demodulator 20. The synchronous demodulator 20 generates an range output 22 
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indicative of the distance between an occupant 24 and a potential impact point within the 
vehicle. The range output 22 is supplied to a signal processor/discriminati n unit 26 for 
storage in a memory 36, and for subsequent use as a basis for selecting the temporally 
most appropriate or "optimal" one of a plurality of predetermined crash discrimination 
strategies. The active infrared position detector 12 and the signal 
processor/discrimination unit 26 receive power from a vehicle battery 32. 

In operation, a narrow infrared beam 34 is directed horizontally at the expected 
position of the driver or passenger. At least a portion of the light energy scattered by the 
occupant's clothing is detected by the optical receiver 18, which is located away from the 
beam axis so that the optical receiver 18 can detect differences in reflected light intensity 
and angle associated with occupant absence, presence and forward motion. The infrared 
beam 34 is distinguished from ambient light of the same wavelength by modulating the 
intensity thereof. A modulation frequency of about 10 Khz or higher provides acceptable 
modulation since a minimum of 10 Khz is well within the range of most existing electronic 
and optical systems, and is fast enough to monitor passenger movement in real time. 
Determination of the instantaneous distance of the occupant from the optical transmitter 
16 is based on the relative intensity of the light scattered by the occupant as measured by a 
detector or detectors, with measurements based on relative intensity or the angle from 
which the light is received. 



20 



25 



Referring to FIGS. 2a~c, the infrared beam 34 is be generated by an infrared light 
emitting diode or semiconductor laser 116 located within the optical transmitter 16 and 
is collimated by a lens 118 to produce an infrared beam 34 that maintains a constant 
diameter of approximately one inch over a distance of 6 to 30 inches. The collimated 
infrared beam 34 is then scattered in all directions by a reflecting surface 120, i.e., the 
vehicle occupant 24 if present. The distance detected by the active infrared position 
detector 12 is based upon triangulation of the received light intensity. The scattered 
collimated infrared beam 34 is focused by an imaging lens 202 to form a relatively small 
spot 204 on a position sensitive detector (PSD) 302. Since the imaging lens aperture 
selects the fraction of the scattered light that reaches the position sensitive detector (PSD) 
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302, the spot 204 on the detector moves when the angle a between the axis of the 
infrared beam 34 and the axis of the reflectively scattered light reaching the detector 
changes, as shown in FIG. 2c. 

The PSD 302, as shown in FIG. 2b, is a distributed photosensitive device for which the 
5 contrast ratio (I^VCIi* I2) in output currents I and O from the inside and outside ends 304 
and 306 respectively provides a linear measure of the spot's vertical position. As shown in 
FIG. 2c, when the angle a varies, the position at which the reflected light is imaged will 
vary across the PSD 302. The ratio of the two current outputs I and O varies as the center of 
light intensity moves across the PSD 302 and, therefore, provides a measure of the angle 
10 a . 

The direction of the reflected rays reaching the PSD 302 is determined by the positions 
of the center of the imaging lens 202 and the point where the reflecting surface 120 cuts 
the infrared beam 34. Thus, the particular angle a at which light is received at the PSD 
302 depends only on the distance x to the reflecting surface. This angle is determined by 

15 tan a =b/x, where x is the distance from the source to the reflecting surface, and b is a 
predetermined lateral separation of the transmitter and the detector. As the angle a varies, 
the relative amounts of current I and O is proportional to the location of the spot 204. 
Signal processing of the detector output currents preferably includes the step of calculating 
the amplitude-independent ratio of currents to correct for reflection variation at the 

20 reflecting surface 120. The ratio of the relative intensities determines the location of the 
spot 204 to provide a good measure of the angle a and, hence, of the distance x to the 
reflecting surface 120. 

Alternately, other types of position sensors may be incorporated in the instant invention 
instead of the active infrared position detector 12 illustrated in Figs. 1 and 2a-c. For 
25 example, any type of non-contacting position sensor that measures position by 
triangulation, imaging, or temporal ranging could be incorporated as the position sensor 
300 of the instant invention. Moreover, the position sens r 300 is not limited by the 
frequency or mode of operation, and accordingly may incorporate optical, electromagnetic 
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or ultrasonic signals. More particularly, the position sens r 300 may comprise a radar or 
an ultrasonic ranging sensor. However, in the preferred mode of operation, the instant 
invention incorporates a position sensor 300 in accordance with the teachings of U.S. 
Patents 5,446,661 or 5,490,069, which are incorporated herein by reference. 

First Embodiment 
Overview : 



10 



15 



Referring to Fig. 3, in accordance with a first embodiment of the instant invention, the 
algorithm incorporated therein uses information obtained from a collimated, modulated, 
optical beam that is reflected from a surface and imaged by a lens on a position sensitive 
detector (PSD) 302. The algorithm also uses information from a seat weight sensor, a 
speed sensor, and a door switch. The position of the imaged light on the PSD 302 allows 
„the algorithm to compute the distance to the reflecting surface. The algorithm bases its 
assessment by keeping track of the distance, patterns of motion and weight over time and 
builds up confidence as time progresses. That is, as the scenario gets 'older', the algorithm 
gets more and more confident about its assessment. The algorithm uses signals from the 
optical sensor which are : 

1 . Inside channel signal ( I ) : It is the signal from the inside end of the PSD 302. The 
signal level gets higher as the distance to the target gets higher. 

20 2. Outside channel signal ( O ) : It is the signal from the outside end of the PSD 302. The 
signal level gets lower as the distance to the target gets higher. 

3. Led signal : This is the value of the driving voltage of the LED 116. 

4. Limit detector signal : This signal from the range limit detector 140 becomes high 
when the target is closer than a threshold. 



25 



Besides the optical signal, the algorithm receives weight information from the weight 
sensor 350 and speed information from the speed sensor 50. The door switch 60 is used 
for instantaneous recognition of new scenarios. 
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The algorithm calculates four measures from the input signals. These measures are 
reset at the onset of a new scenario and therefore the recognition of a new scenario is a key 
factor. Based on these measures the algorithm determines whether the scenario is a Rear 
facing infant seat (RFIS) or an empty seat or a person. The airbag is then enabled or 
5 disabled depending on what the most probable scenario warrants. 

Algorithm : 

The algorithm is best understood with reference to Figs. 4a through 16c. Table 1 lists 
the pertinent parameters of the algorithm, and their associated nominal values, for 
parameters referenced herein. In general, the names of parameters are capitalized, whereas 
l o the names of variables are in written in lower case. 

The overall occupant discrimination algorithm 400 in accordance with the instant 
invention is illustrated in the block diagram of Figs 4a and 4b, which includes the 
illustration of various subprocesses as single blocks therein. More particularly, the 
occupant discrimination algorithm 400 incorporates subprocesses 400, 500, 600, 700, 
15 800, 900, 1000, 1100, 1200, 1300, 1400, 1500, and 1600, which are illustrated in greater 
detail in respective Figs, 5 through 16c. The occupant discrimination algorithm 400 is 
described hereinbelow with respect to the respective subprocesses. 



TABLE 1 



Parameter 


Description 


a,b,c 


Calibration constants of theactive infrared position detector 


ACCDLT 


Dampfactor for calculating accacc which in turn is used to calculate the damped 
acceleration, effmot. Nominal value is 0.9914. 


DCTIM 


If quicktim is greater than this value (seconds) then the scenario flag is set to 1 


EFFTHR 


The threshold to which effmot is compared. If effmot is above this threshold 
then that indicates high acceleration or deceleration and slant and variance 
calculations are halted. The nominal value is 6. This is based on taking 1 6 mph 
change in 6 seconds as the tolerance. 


ELTHR 


If element counter (elcnt) is greater than this threshold (seconds) then p1 and 
n1 are updated. Nominal value is 9 seconds. 


FUZTIM 


When the scenario timer is greater than this value (in seconds), the probability 
measure for delta and variance is initiated. Nominal value is 9 (seconds). 


LMBS 


The airbag is disabled if the probability of disable is greater than this threshold. 
Nominal value is 0.5. 


LSGTTHR 


When the low signal timer (lowsig_timer, which counts for how long the 
algorithm has been in the lowsignal mode) exceeds this value ( in seconds) the 
scenario is reset. Nominal value is 45 (seconds). 
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MAXA 
MAXACC 



MAXRET 



M1ND1S 

RESET_SPElD 
THRESHOLD 



SCTIM 



SDEV CLIP 



Maximum possible distance (cms). Nominal value is 120 cms . 
Maximum acceleration tolerated. The absolute value of 
acceleration/deceleration is clipped to this value. The nominal value is 0 5g 
This value of 0.5g converted to mph per 50 millisecond is 0.559, where it is 

calculated using :0.5g=(500»3600)/(3»1760»30.48»20)=0.559 

The value (in seconds )from which the low_return Jimer (the timer which starts 
counting down to zero once the algorithm comes out of the low signal mode) 
starts counting from. Nominal value is 15 seconds. 

Minimum distance when limit detector is on (cms ). Nominal value is 5 cms 
If the speed is above this threshold (mph) then some of the reset parameters 
are changed so as to make the reset of a scenario more 'difficult'. Nominal 
value is 0 (mph). 



When the scenario timer (scentim) is greater than this value (seconds) the 
scenario flag is set to 1. Nominal value is 45 seconds. 

[Maximum value to which the variance measure is clipped. Nominal value is 
0.03. 



SPDVIB 



incremental factor which linearly increase the overall slant threshold with 
respect to speed, to account for vibration induced through the motion of the 
vehicle. 



WGTDtS 



WGTTHR 
WGTTHR 2 



WIOTH 



Threshold for the moving average of the absolute difference of primary and 
[secondary animacy probabilities above which scenario is reset . This threshold 

is used when the vehicle is in motion, Nominal value is 0.9. 

Weight threshold (lbs) used in the weight reset routine when the vehicle is 
stationary. Nominal value is 34 lbs. 

[Weight threshold (lbs) used in the weight reset routine when the vehicle is in 
motion. Nominal value is 34 lbs. 



The sum of the inside and outside psd channel signals (IO) have to be below 
this threshold to enter the low signal mode. Nominal value is 0.7 volts. 



WLEDTH 
WLSGTM 



The inside psd channel (I) has to be less than this threshold to enter the low 
signal mode. 

IThe led signal has to be greater than this threshold (in counts) to enter the low 

signal mode. The nominal value is 4094 (corresponding to FFE hex) 

tlf the algorithm has been in the low signal mode for longer than this period ( in 
seconds), but the algorithm has NOT been in the low signal mode for longer that 
LSGTTHR, then the algorithm starts using only weight. Nominal value is 30 
seconds. 



WSCENTM 



X1 



The outside psd channel (O) has to be less than this threshold to enter the low 
[signal mode. 



If the algorithm has been in the lowsignal mode and the scenario timer is lower 
than this value ( in seconds) than only the weight measure is used. Nominal 
value is 128 seconds. 
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Distance: 

The signal processor/discrimination unit 26 of the instant invention samples data 
such as at 2 Khz, which samples are integrated to lowpass filter the data and down sample 
to 20 hz. Such a sample rate is high enough to detect even the fastest human motion. 
5 , Referring to Fig. 6 illustrating subprocess 600, the I and O signals are used to calculate the 
distance using the formula: 

X(i) = (c-I+0)/(a-I+b)) 

where a, b, and c are calibration constants specific to each optical sensor unit. If a limit 
detector is on, i.e., if a limit detector signal exceeds a certain threshold denoting a target 
10 within close proximity, then a default distance is set to a minimum value (MINDIS). If the 
calculated distance exceeds a maximum value (MAXA), then the distance is set to this 
maximum value. A 32 point moving average of the distance is calculated (avgdis). 

Slant: 

Referring to Fig. 7 illustrating subprocess 700, the slant processing subprocess is used 
15 for the detection of small amplitude motion. Small amplitude human motions such as 
breathing are very hard to recognize from changes in the distance measures because of the 
noise in the distance signal. To identify these motions the algorithm looks at the variations 
in the total returned signal, I+O = IO. The addition of I and O cancels the thermal PSD 
noise, since that noise is antisymmetric . A 12 point moving average of IO is done to get 
20 sigav. This is used to calculate a 1.5 second slant transform. At 20 hz this is a correlation 
of the sigav sequence with a 31 point sequence { -15... 15}. This slant value is normalized 
by sigav to get the percent variation The absolute value of the normalized factor is then 
multiplied by a distance factor (slant Jactor) to get a result which corresponds to motion. 
The slant Jactor is calculated as follows: 

25 slantjactor = {(avgdis - 50) * SLFAC + 50}/50, 
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Where SLFAC is a number between 0 and 1 . A typical value is 0.7 . An SLFAC of 1 
would correspond to a linear distance multiplier. Because non-PSD noise tends to increase 
somewhat with distance, a SLFAC of less than 1 is preferably used, which gives a weaker 
than linear multiplication and some suppression of the non PSD noise. This final result is 
clipped and finally a 6.4 second moving average is performed. This is the final slant value ( 
slantav). This slantav is compared to a speed dependent threshold (thrcnta), where the 
threshold is calculated as follows : 

thrcnta = min [ SLCLIP, (THRCNT + SPDVIB * speed)] 

where, THRCNT is the base slant threshold (at 0 speed) and SLCLIP is the maximum 
allowable slant value which is usually the value at - 60 miles per hour . SPDVIB is the 
incremental factor which linearly increases the overall slant threshold (thrcnta) from a 
value of THRCNT at zero velocity to SLCLIP at its specified speed. This is to take into 
account the vibration induced through the motion of the vehicle. If slantav is above thrcnta 
it denotes motion thus indicating the presence of a person. If slantav is below the threshold 
that indicates inanimacy, therefore indicating the presence of either a RFIS or an empty 
seat. 



20 



25 



Speed and damped acceleration: 

Referring to Fig. 8 illustrating subprocess 800, at first a 12 point moving average of 
the speed signal is taken. Acceleration is calculated from the filtered speed signal by taking 
the difference between consecutive samples. The absolute value of the acceleration signal 
is taken and clipped to a maximum tolerable acceleration ( MAXACC) to obtain absacc. A 
moving sum of the absolute acceleration (absacc) is calculated (accsum). The window 
length is MACCN. Any older value outside the window is damped and added to the 
moving sum to obtain effmot. The damp factor is ACCDLT. Therefore effmot is : 

» n- MACCN 

effmot(n) = £ absocc(i) + £ (ACCDLT)> absacc(n - MACCN + 1 - /) 
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Effmot is used as an indicator of high acceleration and deceleration. If it exceeds a 
threshold (EFFTHR) then the motion measures are disregarded. 

Big change : 

Referring to Fig, 9 illustrating subprocess 900, when two consecutive I+O signal 
5 samples go through a big change [ the absolute difference between IO(n) and sigav(n-l) is 
greater than 12.5% of sigav(n-l) for two consecutive samples ] and/or the limit detector 
stays on for two consecutive samples then a 'big change ' is said to have taken place. A 
timer (quicktim) starts timing at this point. This timer tells when the last ' big change' 
occurred. 

10 Algorithm measures : 

Referring to Figs. 10 and 11 illustrating subprocesses 1000 and 1100, which 
respectively calculate primary and secondary measures, the algorithm calculates four 
measures from the optical and weight signals. The measures are 

1 . Mean distance( x ) 
is 2. Variance of the distance (a 2 ) 

3. Mean weight (w) 

4. Slant measure ( S ) 

Two sets of measures are calculated. One is called the 'primary' measures which are 
calculated from the beginning of the scenario and is reset at the onset of a new scenario. 
20 The other is the 'secondary' measures which are calculated from the onset of a 'big change' 
and are reset at the onset of the next 'big change'. 

Mean distance(x ) : 

This is the mean of the distance of the scenario calculated from the very beginning of 
25 the scenario,ie. 
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1 n 

x ( n ) = ~ 2]avgdis(i) 



where, n is the total number of samples available since the beginning of the 



scenario. 



Similarly, the secondary mean distance is calculated, where n is the number of samples 
from the last big change 

Variance of the distance (cs 2 ) : 

This is a measure of the large amplitude motion. Once the scenario ' settles down' ( the 
scenario flag is set) and the limit detector is not on and there is no high acceleration or 
deceleration present ( effmot < EFFTHR) a running variance of the distance is calculated. 
The distance is compared to a moving average distance ( X v ) as opposed to the mean 
distance till the present instant. The variance is not updated when effmot > EFFTHR or if 
the limit detector turns on because in these cases the signal level would show variance 
which are not induced by the actual motion of the target, but due to the vehicle 
acceleration/deceleration or the blocking of the optical sensor. The variance is calculated as 



15 a 2 (n) = - £ { avgdis(i) - X v (i)} 2 



n ... 



The variance is clipped to a maximum value (SDEVCLIP). 

Similarly, the secondary variance is calculated, where n is the number of samples from 
the last big change. 

Mean weight of the scenario: 

20 The scenario mean weight is taken from the beginning of the scenario. 

iv(n)= ^2>(i) 



n i-. 
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The secondary average weight is calculated in a similar way, where n is the number of 
samples from the last big change. 

Slant measure (&): 

Two motion counters pi and nl are used to indicate animacy and inanimacy 
5 respectively. If the slant value (slantav) is above a threshold (thrcnta) then pi is 
incremented , otherwise nl is incremented. The counters pi and nl are reset at the onset of 
a new scenario. There are similar counters displ and disnl which are analogous to pi and 
nl respectively. The only difference is that these are reset every time there is a big change. 
The counters pi and nl are not updated till the scenario flag has been set for a certain 
10 period. This is checked by seeing if elcnt is greater than ELTHR. This is to take into 
account the filter delays associated with the slant value. The second condition for not 
updating pi and nl is that if the limit detector is on. The pi and displ counters are not 
updated in the presence of acceleration or deceleration. The final motion measure is 

6= pl/(pl+nl) 

*■ 

15 This is the primary motion measure. The secondary motion measure is : 

5 2 = displ /(displ + disnl) 

This measure has a range of from 0 to 1 . The lower the value, the more it is indicative 
of an inanimate object (RFIS or empty seat). 

Low Signal Routine 1200 : 

20 Referring to Fig. 12 illustrating subproccss 1200, the algorithm enters a low signal 

mode if any one or more of the following three occur : 

1) The inside channel (I) is below a threshold ( WITH) 

2) The outside channel (O) is below a threshold (WOTH) 
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3) The sum of the inside and outside channels (IO) is less than a threshold (WIOTH) 
and the LED signal is above a threshold (WLEDTH) 

This could occur in the event of a very dark target (eg, black velvet). In this case the 
LED pumps up to almost its maximum value (if not its maximum value) and the returned 
5 signal has a low signal to noise ratio. The algorithm freezes off the I and O signals, ie 
disregards the current I and O signals. In other words : 

I(n) = I(n-l) 
0(n) = 0(n-l) 

At the same time a flag (lowsig_flag) is set to denote the presence of the low signal 
1 0 mode. A timer (lo wsigjimer) starts timing at the same time to keep a track of how long the 
system has been in the low signal mode. 

The algorithm constantly monitors to see if the optical sensor is blocked or the signal 
reflected off the target has a low signal to noise ratio. Once that happens, the algorithm 
'recognizes' that and takes alternate routes and continues its assessment using the available 

15 useful signals. If the algorithm has been in the low signal mode for a long time then the 
algorithm reverts to using just weight. Till that happens, the algorithm still relies on its 
measures derived from the optical signals which are nothing but the last valid values saved 
prior to entering the low signal mode. This is done since the likelihood of the patterns 
changing during that time is low. If, however, the scenario is fairly new, then from the very 

20 onset of the low signal mode the algorithm uses just weight. 

Once the algorithm comes out of the low signal mode, it does not start using the 
variance and delta measures right away. A timer (low_return_timer) starts timing out from 
a set value (MAXRET) and once it times out to zero the two motion measures are used. 
This is analogous to the startup condition. Till this timer times out to zero, the low signal 
25 timer continues to count on. If this counter exceeds a threshold (LSGTTHR) indicating that 
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the algorithm has been in the low signal mode for a long enough duration to warrant 
gathering information afresh. The scenario is reset in this case. 

Once the algorithm enters the low signal mode it freezes off the optical signal to the 
value prior to entering the low signal mode and uses for a little while till that is rendered 
5 too old. On entering the low signal mode the algorithm disables the airbag in the event of 
an empty seat even if the requirement is to enable the airbag with an empty seat. This is 
done because the only useful information at this point is weight and a RFIS could weight 
light too and the scenario may come up to be an empty seat and enable the airbag. 

Scenario settling down: 

10 Referring to Fig. 13 illustrating subprocess 1300, once a new scenario is recognized, 

then a flag (scenflg) is cleared. A timer ( scentim) starts timing. Once the ' big change' 
timer (quicktim) counts up to a threshold (DCTIM) then the scenario flag is set. If the big 
change timer hasn't timed out to the threshold because of frequent big changes but the 
scenario timer has timed out to its threshold (SCTIM) then the scenario flag is set. At this 

15 point a counter (elcnt) starts counting. Thus the difference between this counter and 
scenario timer tells how long it took for the scenario to settle down. The scenario at this 
stage is said to have 'settled down'. This is done in order to disregard motion that may 
occur at the beginning of the scenario. For example, a RFIS may be put in the seat and the 
scenario change is recognized at that point. There will probably be motion as the RFIS is 

20 being positioned properly and then belted up. If this initial motion is considered then the 
motion measure for this scenario would indicate high motion and as the algorithm 
remembers what it sees for a long time, this high motion would have an effect for a long 
time. By letting the scenario to 'settle down', the algorithm disregards all this initial 
motion. 
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Probability calculation 

Referring to Figs. 14a, 14b, and 14c illustrating subprocess 1400, the scenario sample 
space { Sj } is divided into the following three categories: 

(SO) Rear Facing Infant Seat (RFIS) 



(SI) 
(S2) 



Empty Seat (ES) 

« 

Person (PSN) 



10 



All these scenarios are assigned a priori probabilities, ie, the a priori probabilities 
Pr(RFIS) 
Pr(ES) 
" Pr(PSN) 



are 



We have probability density functions defined for the above four measures given these 
three scenario : 



WMJS.) 

where = {x , a 2 , w,5) 
15 and S, = { RFIS, ES, PSN} 

All these scenarios are assigned a priori probabilities Pr(Si). We have 
probability distributions defined for the above four measures given these three 
scenario [UM^ J. So that makes a total of 12 distributions. As the algorithm 
calculates the four measures, probability of a particular measure given a particular 
10 scenario [Pr(Mj|Si) ] is calculated using the available distributions. The probability 
of the four measures, given a particular scenario, is calculated as follows : 
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Pr(M 0 ,M„M 2 ,M J |S i ) = flPk(M|Si)Wji + Z ZPr(M|SOPr(Mc|Si)Pr(S i ){C jkKj } 

j*k, k>j 

where, W jr and C jkK are the autocorrelation and crosscorrelation terms respectively. The 
correlation matrix (l { ) for each scenario for the five measures is : 



10 



li = 







0)2 Ki 


0)3 K, 




w 






^20^ 












^32 Ki 





(i = 0,1,2) 



It is reasonable to assume that the measures are independent. For example, the weight 
and the average distance obviously are independent. The motion measures are also 
independent of the range and weight since they are invariant to these measures. It is, 
therefore, assumed that there is very little correlation between the measures and the cross- 
correlation terms are taken to be zero. The auto-correlation terms being one, the above 
equation then simplifies down to : 



Pr(Mo,M 1 ,M 2l M s |S l ) = FI Pr ( M jl Si ) 

j=0 

Therefore, the probability of the four measures given a particular scenario is : 

Pr { ( x , a 2 , w , 5) | RFIS} = Pr( 3c |RFIS) Pr^RFIS) Pr( w |RFIS) Pr(8|RFIS) 
15 Pr { ( x , a 2 , w , 8) | ES } = Pr( x |ES) Pr(o : |ES) Pr( w |ES) Pr(8|ES) 

Pr { ( x , a 2 , w , 5) | PSN} = Pr(3c |PSN) Pr(a'|PSN) Pr( w |PSN) Pr(8|PSN) 

Then using Baye's rule, the probability of a scenario given the four measures is 
calculated : 



Pr(M o,M .,M 2,M 3|S.)Pr(Si) 
Pr(S i |M 0 ,M 1 ,M 2 ,M 3 ) -— 2 - 



£ Pr(Mo,M i,M 2,M j|Sk)Pr(Sk) 

k - 0 



(4) 



20 



for i = 0,1,2 
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ie, 



Pr{RFIS|(*,cr',w,5)} = 

= Pr{(x,oJ , w,8 )| RFIS) Pr(RFIS) 

Pr{(x,a^w,8)|RFIS}Pr(RFIS) + Pr{(x,a^w,5)|ES}Pr(ES) + P^{(x > a^w,5)|PSN}Pr(PSN) 

(5) 

Pr{ES|(x,a 2 ,iv,S)} = 

_ Pr{(x,a : ,w,6)|ES}Pr(ES) ^ 

Pr{(x '° ; .w,5)|RFIS}Pr(RFlS) + Pr{(x,o 3 , w,8 )| ES> Pr(ES) + Pr{(x,a 2 , w",6)| PSN} Pr(PSN )' 

(6) 

Pr(PSN|(x,a 2 ,vi;,8)} = 

_ ■ Pr {(x,o 2 , w,8)| PSN}Pr(PSN) 

Pr{(x,a^w,8)|RFIS}Pr(RFlS) + Pr{(x,a^w > 5)|ES}Pr(ES) + Pr{(x,a^w,5)|PSN}Pr(PSN) 

(7) 

So we get three probabilities, ie, probabilities of the three scenarios, given the four 
measures, and these tell us which scenario is the most probable. The airbag is then enabled 
or disabled depending on what the most probable scenario warrants. It may be added here 
that a Hag is used to enable or disable with an empty seat. If it is desired to disable the 
airbag for an empty seat, then the sum of the probability of RFIS and empty seat 
determines whether the airbag should be disabled. This probability is compared with a 
threshold and if it is above this threshold then the airbag is disabled. 

The algorithm may be made to run without using one or more measures. For example, 
it may be desired to see how well the algorithm performs in its detection capability with 
only the optical sensor. In this case a flag is used to disable the use of the weight measure 
in the algorithm. In that case, the probabilities of weight given the individual scenarios are 



25 set to 1 , ie 



Pr( vv | RFIS) = Pr( w |ES) = Pr( w |PSN) = 1 



Equations (1), (2) and (3) therefore reduce to : 
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Pr { {x , a 2 , 5) | RFIS} = Pr(3c |RFIS) Pr(cr|RFIS) Pr(5|RFIS) 



(8) 



Pr { (x , a 2 , 8) | ES} = Pr( 3c |ES) Pr(o~'|ES) Pr(5|ES) 



(9) 



Pr{(3c, a 2 , 6)|PSN} = Pr( x |PSN) Pr(a'|PSN) Pr(8|PSN) 



(10) 



The probabilities of the scenarios thus reduce to : 



Pr{RFIS|(x,a 2 ,8)} = 



15 



Pr{(x,a 2 , 5)| RFIS} Pr(RFIS) 



Pr {(x,a 2 ,5 )| RFIS} Pr(RFlS) + Pr{(x,a 2 ,8 )| ES} Pr(ES) + Pr {(x,o 2 ,8 )| PSN } Pr(PSN) 



io Pr{ES|(3c,a 2 ,8)} = 



Pr{(x,a 2 ,S)|ES}Pr(ES) 



Pr{(x,o 2 ,8)|RFlS}Pr(RFIS) + Pr{(x,o 2 ,6)|ES}Pr(ES) + Pr{(x,a 2 ,8)|PSN}Pr(PSN) 



Pr(PSN|(x ,a 2 ,8)} = 



Pr{(x,o 2 ,8)|PSN}Pr(PSN) 



Pr{(x,o 2 , 8)| RFIS} Pr(RFIS) + Pr{(x,a\S)|ES}Pr(ES) + Pr{(x,o 2 ,8)|PSN}Pr(PSN) 



(11) 



(12) 



(13) 



The probability based on distance and weight are calculated from the very beginning of 
20 the scenario. However, the variance and motion measures are not used right away for 
calculating probabilities. The scenario timer (scentim) has to be above a threshold 
(FUZTIM) to initiate the calculation of the probability based on variance and motion 
measures. Till then, 
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Pr(a'IRFIS) = Pr(<r|ES) = Pr(o~'|PSN) = Pr(5|RFIS) = Pr(8|ES) = Pr(6|PSN) = 1 
S condary Probability c Iculation 

Referring to Figs. 15a, 15b, and 15c illustrating subprocess 1500, the secondary 
probabilities are calculated in exactly the same manner-the only difference is that the 
measures used are the corresponding secondary measures. The other difference is that all 
four probabilities are calculated from the very beginning of the period between two 
consecutive big changes-there is no wait period for the initiation of the variance and 
motion measures because they are on going even through the occurrence of a big change. 

Recognition of a New Scenario : 

Referring to Figs. 16a, 16b, and 16c illustrating subprocess 1600, a new scenario can 
be detected using any one or more of the following five conditions : 

1 ) Door switch :In most, if not all, cases the passenger side door will be opened and/or 
closed while changing a scenario. Therefore the door switch is used as the primary 
indicator of a change in a scenario. 

2) Weight reset : The change in a scenario will also accompany a change in weight 
and/or distance. Therefore if the absolute change in weight and distance exceeds a 
threshold and this change persists for a while then the scenario is reset. The weight 
and distance are compared to moving average values of weight and distance so that 
an instantaneous change in either does not initiate the change in a scenario. 

3) Resets using probabilities : This is where the secondary probabilities are used. If 
there has been a change in the scenario, then there will be a persistent difference in 
the primary and secondary probabilities of a particular scenario. This is because, the 
primary probabilities are calculated from the primary measures which have a long 
time history and are less prone to change in short duration change in pattern of the 
measures. The secondary probabilities, on the other hand, are measured from the 
secondary measures, which have a shorter time history and therefore track real 
changes in the measures quicker. If there has been an actual change in the scenario, 
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then the secondary measure would persist to be indicative of the new scenario, 
whereas the primary would still tend to 'hold' on to its older values. There would 
be a persistent difference in the primary and secondary probabilities in this case. If 
this change persists for a certain duration of time then the scenario is reset. The way 
5 this is done is that the absolute difference of the primary and secondary 

probabilities are taken and are integrated. If this integrated value exceeds a 
threshold then the scenario is reset. Resets AA, BB and CC in the algorithm are 
these resets in the cases of RFIS, empty seat and people respectively. 

The same set of reset conditions are checked with different thresholds in the event of 
10 the vehicle being in motion. The reason behind doing this is that it is much less likely for 
the scenario to change while the vehicle is in motion than that when the vehicle is 
stationary. Therefore if the vehicle is stationary the resets are made quicker and when the 
vehicle is in motion the resets are made harder. 

Once a new scenario is detected all the measures, counters are reset. 
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Second Embodiment 

A second embodiment of an occupant discrimination algorithm in accordance with the 
instant invention is illustrated in the block diagram of Fig. 17. Table 2 lists the pertinent 
parameters of the algorithm, and their associated nominal values, for parameters referenced 
herein. 



TABLE 2 



ANIMULT 


0.156 


ANITHR 


0.55 


BSMULT 


0.01 


DECSLW 


0.984 


DISBS 


40 cm 


DISSUN 


16 cm 


FIRETIM 


16 cm 


FSTBEG 


4 sec 


FSTEND 


8 sec 


NITHR 


80 


PRBS 


0.5 


PRNWS 


0.55 


PRSLIN 


0.1 


QCKTO 


45 sec 


QTBS 


1 20 sec 


QTNWS 


60 sec 


QUICK 


3.2 sec 


RECDEC 


0.938 


SCENBS 


1 20 sec 


SEATBCK 


65 cm 


SEATCHG 


8 cm 


THINDIS 


20 cm 


THINSL 


0.3 


THRSLW 


240 



The following calculations are made on an ongoing basis: 



10 1. Distance: 

The distance calculations in subprocess 600* are the same as subprocess 600 of the first 
embodiment. 
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2. Slant: 

This is used for small amplitude motion detection. Small amplitude human motions 
such as breathing are very hard to recognize from changes in the distance measures because 
of noise in the distance signal. To identify these motions, the signal 

5 processor/discrimination unit 26 examines the variations in the total returned signal, 
I+Q=IO. Because the total signal varies inversely with the square of the distance, using the 
total signal increases signal to noise ratio by a factor of 2. The addition of I and O also 
cancels any thermal PSD noise, since that noise is antisymmetric (ie, positive noise 
excursions in I = negative noise excursions in O ). This accounts for about 50% of the 

10 noise, thus producing an additional doubling of the signal to noise ratio. Furthermore, 
many human motions result in changes in the reflectance, such as created by changes in the 
reflection angle or color of the reflecting surface. This can further enhance the SNR. 

The 1 and O signals are added to produce 10. A 6 point moving average of IO is done 
to determine sigav. This is then used to calculate a 1.5 second slant transform. At 20 hz, 
15 this is a correlation of the sigav sequence with a 31 point sequence {-15... 15}. The slant 
value is normalized by sigav to get the percent variation. The absolute value of the 
normalized factor is then multiplied by a distance factor (slantjactor) to get a result which 
corresponds to motion. The slantjactor is calculated as follows: 

slantjactor = {(avgdis - 50) * SLFAC + 50}/50, 

20 where SLFAC is a number between 0 and 1. A typical value is 0.7. An SLFAC of 1 

would correspond to a linear distance multiplier. Because non-PSD noise tends to increase 
somewhat with distance, an SLFAC of less than 1 is preferably used, which gives a weaker 
than linear multiplication and some suppression of the non-PSD noise. This produces the 
final slant value (slantav). 



* » 
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Once certain conditions are satisfied as explained hereinbelow, the slant value is 
compared to a threshold for calculating the slant measure ( ). If slantav is greater than this 
threshold and there has not been high acceleration, then a counter pi is incremented, 
otherwise a counter nl is incremented. 



5 3. Speed and damped acceleration: 

At first, a 12 point moving average of the speed signal is taken. Acceleration is 
calculated from the filtered speed signal by taking the difference between consecutive 
samples. The absolute value of the acceleration signal is taken and clipped to a maximum 
tolerable acceleration (MAXACC) to obtain absacc. A moving sum. of the absolute 
10 acceleration (absacc) is calculated (accsum). The window length is MACCN. Any older 
value outside the window is damped and added to the moving sum to obtain effmot. The 
damp factor is ACCDLT. Thus, effmot is: 

effmot(m) = 1/K * [accsum(m) + accacc(m)], 

where accacc(m) = ACCDLT * [accacc(m-l) + absacc(m-MACCN)]. 

1 5 Conditions: 

The algorithm constantly checks for the following conditions: 

1 . Big change: 

When two consecutive I+O signal samples go through a big change, e.g., the absolute 
difference between IO(n) and sigav(n-l) is greater than 12.5% of sigav(n-l) for two 

20 consecutive samples, and/or the turning on or off of a limit detector persists for two 
consecutive samples, then a "big change" is determined to have taken place. Two timers, 
i.e., quicktim (0.05 second increments) and dischtim (6.4 second increments), are started. 
These timers tell when the last "big change" occurred. At the same time, a set of motion 
counters sidnl, displ, disn2 and disp2 are updated. These motion counters follow the same 

25 criteria as those of updating nl , pi , n2, p2. 
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When the next "big change" is encountered, the timers and the counters are reset. 
Another counter (bccnt) is incremented every time there is a "big change," thereby tracking 
the number of times a "big change" has been detected in a scenario. This counter is reset at 
the onset of a new scenario. The 32 point moving averaged distance (avgdis) prior to the 
5 big change is saved. 

2. Scenario settling down: 

Once a new scenario is recognized, then a flag (scenario flag) is cleared. A timer 
(scenario timer) starts timing. Once the "big change" timer (quicktim) counts up to a 
threshold (DCTIM), the scenario flag is set. If the big change timer has not timed out to 
10 the threshold because of frequent big changes, but the scenario timer has timed out to its 
threshold (SCTIM), i.e., the new scenario was recognized SCTIM seconds ago but the last 
big change was noted less than DCTIM seconds ago, then a scenario flag is set. The 
scenario at this stage is said to have "settled down." 



Normal Operation: 

1 5 Normal Operation is outlined in the flow chart of Fig. 2. As noted above, the five 

measures calculated during the discrimination process of the instant invention are as 
follows: 

(M 0 ) Mean distance of a scenario (avgda): From the inside (I) and outside (O) 

signals the distance (X) is calculated. The mean of the distance of the scenario 
20 is calculated from the very beginning of the scenario, i.e.. 

1 ■ 

avgda(n) = — ffcAXO 

n ui 

(M,) Maximum distance (X max ): The maximum distance since the beginning of the 
scenario is recorded; 
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(M2) Variance of the distance (a 1 ): This is a measure of the large amplitude motion. 
Once the scenario "settles down," a limit detector does not indicate a close 
proximity (i.e. "turns on"), and there is no high acceleration or deceleration 
present (effmot < EFFTHR), a running variance of the distance is calculated. 
The distance is compared to a moving average distance (Xv) as opposed to the 
mean distance till the present instant. The noise portion (avgda2 * NLP))2 is 
subtracted off to get the actual variance. The variance is not updated when 
effmot> EFFTHR or if the limit detector turns on, because in these cases the 
signal level would show variance which is not induced by the actual motion of 
the target, but due to the vehicle acceleration/deceleration or the blocking of the 
optical sensor. The variance is calculated as: 

2 (n) = max[l/n lon {X(i) - X(i)} 2 - {avgda(n) J *NLP} 2 ,0]. 

(M 3 ) Mean weight of the scenario: A running average of the scenario mean weight 
is taken from the beginning of the scenario: 

15 W(n)=l/nI i . llon w(i) 

(M 4 ) Slant measure ( ): =pl/(pl+nl); 

Once a new scenario is detected, all measures are initialized. The scenario flag is 
cleared. From then on the scenario mean distance, the maximum distance and the scenario 
mean weight are routinely updated. These three measures are used at this point to calculate 

20 the probability of a scenario. The variance and the slant measures are not used at this point 
because at the beginning of a new scenario, there may be motion (for inanimate objects) 
which should be disregarded. At this point, the marginal probabilities for slant and 
variance, ie., Pr( S f ) and Pr( 1 S t ), are set to one. When the scenario "settles down," i.e., 
the last big change took place more than a predetermined period of time ago (DCTIM) or 

25 the scenario is a predetermined period old (SCTIM), whichever of these two takes place 
first, the scenario flag is set and the variance calculation begins. Once the scenario timer 
exceeds a threshold (FUZTIM), the probability based on variance is calculated, i.e., 
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Pr( 2 Sj). The probability of variance given a scenario is calculated using the following 
formula: 

Pr( 2 Si)<-Pr( 2 S,)*|i + (l-»i)»Pi< 2 S 6 ); 

where n = min [ 1, (pi + nl)/PNDIV], and Pr( 2 S 6 ) is the probability of variance 
5 given a scenario where a person is reading a newspaper. This is used as the most universal 
distribution because it is the flattest. 

Once the scenario settles down, an element counter (elcnt) starts counting, and when 
the count exceeds a threshold (ELTHR), the slant measure updating is initiated and the 
probability based on the slant measure is calculated. 

At this point the slant transform (slantav) is compared to a speed dependent threshold 
(thrcnta), where the threshold is calculated as follows: 

thrcnta = min[SLPLIP, (THRCNT + SPDVIB * speed)] 

where, THRCNT is the base slant threshold (at O speed) and SLCLIP is the maximum 
allowable slant value which is usually the value at approximately 60 miles per hour. 
SPDVIB is the incremental factor which linearly increase the overall slant threshold 
(THRCNTA) from a value of THRCNT at zero velocity to SLCLIP at its specified speed. 
This accounts for vibration induced through the motion of the vehicle. 

The motion counters pi and displ are incremented if the slant transform (slantav) is 
greater than THRCNTA and at the same time there is no detection of high acceleration or 
20 deceleration (effmot < EFFTHR). Following the same criterion, two accumulators (recani 
and slowani) are incremented using damp factors, i.e., rcani = recani*RECDEC + 1, and 
slowani = slowani* DEC SLW + 1, where RECDEC and DECSLW are damp factors. These 
accumulators record the motion history of the scenario and gradually damp out older 
information. They are also used in detecting the onset of new scenarios. 
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The counters nl and disnl are incremented if the slant transform is less than 
THRCNTA regardless of the acceleration or deceleration of the vehicle. Following the 
same criterion, recani and slowani are decremented using damping recami = 
recani*RECDEC - 1, and slowani = slowani'DECSLW - 1. The slant measure (symbol) 
5 equals pl/(pl+nl). The slant probability is also calculated like the variance probability as 
follows: 

Pr( S f )<-Pr( S^u + O-jO'Prf S ,); 

where » = min[ 1, (pl+nl)/PNDIV], and Pr( S { ) is the probability of slant measure 
given the scenario is a front facing infant seat. Since the distribution is uniform, this 
10 probability is a constant. 

Signal processor/discrimination unit 26 routinely monitors to determine if the optical 
sensor is blocked or the signal reflected off the target has a low signal to noise ratio. Once 
that happens, signal processor/discrimination unit 26 operates in one of a plurality of a 
fall back modes to continue assessing the other factors using the available useful signals. 
15 These fall back modes are as follows: 

1. Blocked sensor: 

When active infrared position detector 12 is blocked, signal 
processor/discrimination unit 26 continues situation assessment with the remaining 
available signals. If active infrared position detector 12 is blocked from the very 
20 beginning of the scenario, then the assessment is based on the slant measure, variance and 
weight. The marginal probabilities based on distance and maximum distance are set to 
one, thus de-emphasizing them. If active infrared position detector 12 becomes blocked 
sometime within the scenario, then the maximum distance is also used. In this situation, 
only the marginal probability of the distance is set to one. 
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2. Low Signal; 

In the event of a low returned signal as may happen with the target object having a 
highly non reflecting surface, signal processor/discrimination unit 26 becomes "blinded" 
as the signal to noise ratio becomes low. The scenario mean distance, the scenario 

5 maximum distance, the scenario variance and the slant measures are saved. After the onset 
of the low signal mode, for a certain period these four saved measures are used along with 
the scenario mean weight. After that period is over, the saved information is deemed too 
old and signal processor/discrimination unit 26 continues situation assessment using 
only weight. In this situation, the airbag is disabled with an empty seat even if under 

10 normal circumstances the airbag is required to be enabled with an empty seat. This is due 
to the fact that a rfis and an empty seat would both be in the same weight range. 

More specifically, if the LED value exceeds a threshold (WLEDTH) and the IO signal 
falls below a threshold (WIOTH), then the returned optical signal has a low signal to noise 
ratio. Let us assume this is the first time the low signal mode is being entered. A flag 
15 (lowsignal_flag) is set to denote that signal processor/discrimination unit 26 is in the low 
signal mode. 

A timer (lowsig_timer) starts timing to tell how long signal processor/discrimination 
unit 26 operates in the low signal mode. The information prior to entering the low signal 
mode is saved. If the scenario is fairly new, then the saved information is deemed to not be 

20 "mature" enough, i.e., substantial time history has not been built up. This is determined if 
the scenario timer is lower than a threshold (WSECNTM). In such a situation, only the 
scenario average weight is used to calculate the overall probability. The other marginal 
probabilities are set to one. Also, the same is done if signal processor/discrimination 
unit 26 has been in the low signal mode for a long time, in which case all saved 

25 information is deemed to be too old. This is determined by lowsig_timer exceeding a 
threshold (WLSGTM). 
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When the low signal period is over, lowsignal_flag is cleared and a timer 
(lo\y_retum_timer) starts timing out from a maximum value (MAXRET). This timer 
permits the slant information to build up again, i.e., the slant array is filled with pertinent 
samples. A flag (low_ret_flag) is set to denote the "recovery period," while the 
5 lowsig timer continues to count. 

Signal processor/discrimination unit 26 continues its operation by comparing the 
lowsig_timer against a threshold (LSGTTHR). Once this threshold is exceeded, the 
scenario is reset. If not, then signal processor/discrimination unit 26 continues its 
operation using all the measures. During the "recovery period," the low_retum_timer is 
0 decremented and if it times out before the lowsig_timer exceeds LSCTTHR, then the 
following checks are made to make sure that the scenario has not changed. At this time 
low_ret_flag is cleared to denote the end of the "recovery period". A new scenario is 
considered to have taken place if any one of the following occurs: 



15 



20 



25 



1. the scenario is animate (Probability of animacy {Pr(Animate)} is greater than a 
threshold {ANITHR}) and one of the motion history measures (Aniquick) denotes 
inanimacy, i.e., Aniquick is less than a threshold (sigavt * ANIMULT), where sigavt 
is the saved sigav value prior to entering the low signal mode, and ANIMULT is a 
multiplying factor); 

2. The scenario is inanimate [Pr( Animate) < ANITHR] and there has been a big change 
in the distance since the onset of the low signal mode, i.e., the absolute difference 
between the saved distance prior to entering the low signal mode and the 32 point 
moving average distance (avgdis) is greater than a threshold (DISSUN); 

3. The scenario is inanimate [Pr(Animate) < ANITHR] and there has not been a big 
change in the distance since the onset of the low signal mode, i.e., the absolute 
difference between the saved distance prior to entering the low signal mode and the 
32 point moving average distance (avgdis) is less than a threshold (DISSUN), and 
the motion history measure (Aniquick) indicates motion, i.e., Aniquick > sigavt 
♦ANIMULT. 
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Finally, as discussed above, an important factor is to recognize a new scenario and 
restart all calculations. A key trigger indicating a new scenario is operation of the door 
switch, because in most cases the passenger side door will be opened and closed. 
However, in the low probability event of the passenger side door not being used, then 
signal processor/discrimination unit 26 uses the following reset paths. As assessment 
process builds up confidence over time, a determination that a new scenario has occurred is 
not made simply because there has been a sudden change in the pattern of the signals. For 
example, if signal processor/discrimination unit 26 has detected a baby seat for ten 
minutes, and the driver shifts the baby seat momentarily, this is not indicative of a new 
scenario. 

However, if the baby seat were actually being taken out from the seat without opening 
or closing a door, say, in the case of a convertible, then this change in the pattern of motion 
would persist, and signal processor/discrimination unit 26 would eventually decide that a 
new scenario has occurred. 

The time it takes to recognize new scenarios varies depending on which reset path 
initiated the new scenario. Each of these paths are governed by multiple checks which arc 
in turn controlled by adjustable parameters or thresholds. The reset paths are as follows: 

Reset A: 

The scenario is reset if the following conditions are satisfied: 

1 . A "big change" 1 has occurred; 

2. The seatback flag (stbkchg) is cleared; 

3. A certain time has elapsed since the last "big change" (quicktim > QUICK); 

4. The distance is greater than the seatback distance (avgdis > SEATBCK); 

5. Some motion history has built up (recani > NITHR); 

6. The scenario has existed for a while (scenario timer > SCENLM); and 

7. The probabilities show inanimacy (Pr(person) <= PERLM). 

Reset B: 
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A new scenario is recognized if: 

1 . A big change has occurred; 

2. The seatback flag (stbkchg) is cleared; 

3. A certain time has elapsed since the last big change (quicktim > QUICK); 
5 4. The distance is greater than the seatback (avgdis > SEATBCK); and 

5. There has not been much motion built up (recani < NITHR). 

ResetB f : 

A new scenario is recognized if: 

1 . A certain time has elapsed since the last big change (quicktim = QCKTO); 
10 2. The distance is greater than the seatback (avgdis > SEATBCK); 

3. There has not been much motion recently (recabu < NITHR); and 

4. The seatback flag (stbkchg) is set. 

Reset C: 

A new scenario is recognized if: 

15 1 . A certain time has elapsed since the beginning of the scenario (scenario timer > 

SCENBS); 

2. The distance is less than the average baby seat distance in the rear facing mode 
(avgdis <=DISBS); 

3. The probability of rear facing infant seat is low (Pr(rfis)<PRBS); 

20 4. Inanimacy has been noted since the last big change (disnl>BSMULT*displ); 

5. A certain time has elapsed since the last big change (quicktim>QTBS); and 

6. The current scenario is not a person reading a newspaper (Pr(PRN) < PRNWS). 

If conditions 4 or 5 are not satisfied but quicktim > QTNWS, then the scenario is reset. 
This reset path is a "fail safe" reset to correct a situation in which a rear facing infant seat 
25 has been incorrectly classified as something else. ~ 



BNSOOCIO: <WO__93*903i Al I > 



WO 98/49031 



38 



PCT/US98/08429 



Reset D: 

A new scenario is recognized if: 

1 . A big change has occurred; and 

2. The seatback flag is set (indicating a recent sighting of the seat back). 
5 This reset path is to recognize a change from an empty seat. 

Reset E: 

When a certain time has elapsed since the last big change (quicktim=FSTBEG), the 
current filtered 1=0 value (sigav) is stored to be used as a reference value (sigavref) to be 
used in this pass. The motion indication "aniquick" is cleared. From this instant on until 
10 quicktim is equal to FSTEND, aniquick is incremented with the absolute difference of the 
current sigav value and the reference sigav value. At the end of this period, a new scenario 
is recognized if: 

1 . The distance is greater than the seatback (avgdis > SEATBCK); 

2. The filtered distance (avgdis) is greater than the distance prior to the big change by a 
15 threshold (SEATCHG); 

3. Aniquick is less than a threshold (ANIMULT*sigavref); and 

4. The scenario is at least 12.8 seconds old. 

In this case, the seatback flag (stbkchg) is set. 
Reset F: 

20 A new scenario is recognized if : 

1 . The scenario timer is greater than a threshold (FIRETIM); 

2. If the slant measure ( ) denotes inanimacy ( <THINSL); and 

3. There has been a substantial change in the distance ( avgdis - avgda >= 
THINDIS). 
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Reset G: 



A new scenario is recognized if: 



1 . The scenario timer is greater than a threshold (FIRETIM); 

2. The motion accumulator slowani denotes animacy (slowani >=THRSLW); and 

5 3 . The current state of the probabilities denote inanimacy (Pr(inanimacy) > PRSLIN). 

Weight reset: 

If there is a big change in weight as well as the distance persisting for a substantial 
period of time, then the scenario is reset. Here the scenario mean weight and distance are 
compared to moving averaged weight and distance so as not to reset with instantaneous 
10 changes in the weight and distance. The change in weight and distance are compared to 
their respective thresholds (WGTTHR and WGTDIS respectively). 

4 

A decision on whether to disable the airbag(s) is made at the end of every cycle is 
based on the calculated probabilities. The probability of scenarios which warrant the 
disabling of the airbag are added. If this sum exceeds a threshold (LMBS), then the airbag 
15 is disabled. Otherwise the airbag remains enabled. Depending on whether enabling or 
disabling an airbag is required where there is an empty seat, that probability is included or 
excluded in the above sum. Even if the requirement is to enable with an empty seat, the 
airbag is disabled with an empty seat if the system is in the lowsignal mode. 



While specific embodiments have been described in detail in the foregoing detailed 
description and illustrated in the accompanying drawings, those with ordinary skill in the 
art will appreciate that various modifications and alternatives to those details could be 
developed in light of the overall teachings of the disclosure. Accordingly, the particular 
arrangements disclosed are meant to be illustrative only and not limiting as to the scope of 
the invention, which is to be given the full breadth of the appended claims and any and all 
25 equivalents thereof. 
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WE CLAIM: 

1 . A system for discriminating the condition of seat occupancy by an object on a vehicle 
seat and for controlling the activation of a safety restraint system responsive thereto, 
comprising: 

a) a position sensor for sensing the distance between a fixed structure within the 

5 vehicle and the object on the vehicle seat, whereby said position sensor generates at 

least one distinct component of a signal space; 

b) a seat weight sensor for generating a weight signal responsive to the weight of an 
object on the seat, whereby said weight signal constitutes a distinct component of 
said signal space; and 

10 c) A processor operatively coupled to said ranging sensor and said seat weight sensor 

for calculating the likelihood of each of a plurality of seat occupancy scenarios 
from said distinct components of said signal space, and for controlling the 
activation of the safety restraint system responsive thereto. 
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